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ABSTRACT 

In  this  paper,  a  detailed  analysis  is  carried  out  to  study  the  flow  and  heat  transfer  characteristics  of  visco¬ 
elastic  fluid  over  an  exponentially  stretching  sheet  in  a  porous  medium  with  viscous  dissipation.  The  governing  partial 
differential  equations  are  converted  into  non-linear  ordinary  differential  equations  using  suitable  similarity 
transformation  and  then  non-linear  differential  equations  have  been  solved  by  using  Homotopy  Analysis  Method 
(HAM),  which  provides  the  best  possible  way  to  control  and  adjust  the  convergence  region  using  a  non-zero  auxiliary 
parameter  h.  The  effect  of  physical  parameters  on  velocity  and  temperature  is  discussed  through  graphs  explicitly. 
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INTRODUCTION 

In  recent  years,  the  boundary  layer  flow  of  non  -  Newtonian  fluids  over  a  stretching  sheet  plays  an 
important  role  in  applications  in  industries  and  in  many  engineering  problems  such  as  metallurgical  process,  a 
polymer  extrusion  process  involving  cooling  of  a  molten  liquid  being  stretched  into  a  cooling  system. 
The  pioneer  work  of  Sakiadis  [1]  analysed  the  boundary  layer  flow  over  a  continuously  moving  solid  surfaces  of 
infinite  length  and  moving  with  constant  speed.  Crane  [2]  who  extended  the  work  of  Sakiadis  [1]  on  over 
stretching  sheet  such  that  the  velocity  of  the  surface  to  vary  linearly  with  the  distance  from  the  slit.  Later,  many 
researchers  laid  their  attention  on  the  study  of  boundary  layer  flow  of  Newtonian  and  non  -  Newtonian  fluids  over 
different  stretching  sheets,  i.  e.  linear,  quadratic,  power  law  and  non-isothermal  stretching  sheet  [3-6]. 

The  study  of  heat  transfer  of  non-Newtonian  fluid  flows  through  a  porous  medium  has  gained 
considerably  in  the  industrial  applications  of  petroleum  extraction,  filtration  process  and  separation  process  in 
chemical  industries  and  many  others.  Initially,  Cheng,  p  and  Minkowycz  [7]  analyzed  free  thermal  convection 
from  a  vertical  plate  embedded  in  a  saturated  porous  medium,  considering  wall  temperature  as  a  power-law 
function  of  the  distance  along  the  plate.  Consequently,  Subhash  Abel  and  Veena  [8]  Examined  heat  transfer  of  a 
visco-elastic  fluid  flow  through  a  porous  medium  adjacent  to  a  stretching  sheet  including  the  effects  of  frictional 
heating  and  internal  heat  generation  or  absorption.  Subsequently,  diffusion  of  chemical  reactive  species  of  a  non- 
Newtonian  fluid  was  investigated  by  Prasad  et  al.  [9]  through  the  porous  medium  over  a  stretching  sheet. 
Krishnambal  and  Anuradha  [10]  examined  radiation  effects  on  heat  transfer  of  visco-elastic  fluid  flow  through 
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porous  medium  past  a  stretching  sheet.  M.  Subhas  Abel,  S.  K.  Khan  and  K.  V.  Prasad  [11]  studied  combined  buoyancy 
effects  to  flow,  heat  and  mass  transfer  in  a  visco-elastic  fluid  through  a  porous  medium  over  a  stretching  sheet. 

The  problem  of  viscous  dissipation.  Joule  heating  and  heat  source/sink  on  non-Darcy  MHD  natural  convection 
flow  over  a  flux  permeable  sphere  in  a  porous  medium  is  numerically  analyzed  by  Yih  [12].  The  work  of  Sonth  et  al.  [13] 
analyzed  the  effect  of  the  viscous  dissipation  term  along  with  temperature  dependent  heat  source/sink  on  momentum,  heat 
and  mass  transfer  in  a  visco-elastic  fluid  flow  over  an  accelerating  surface.  Later,  Chen  [14]  examined  the  effect  of 
combined  heat  and  mass  transfer  on  MHD  free  convection  from  a  vertical  surface  byohmic  heating  and  viscous  dissipation. 
Recently,  Mahentesh  M.  Nandeppanavar  [15]  considered  the  flow  and  heat  transfer  characteristics  of  a  visco  -  elastic  fluid 
through  a  porous  medium  over  an  impermeable  stretching  sheet  with  viscous  dissipation. 

The  solutions  given  for  stretching  sheet  in  the  above  literature  were  not  unique  and  derived  another  closed  form  of 
solutions.  Consequently,  Troy  et  al.[16]  investigated  momentum  flow  of  visco-elastic  fluid  containing  exponential  terms  in 
boundary  conditions.  The  results  obtained  were  physically  realistic  produces  a  boundary  layer  only  slightly  altered  in  its 
dimensions  from  the  viscous  one.  In  view  of  the  above  discussions  on  boundary  conditions,  we  present  in  the  next  section 
the  physically  realistic  sequential  similarity  solutions  of  visco-elastic  boundary  layer.  In  view  of  this  point, 
Elbashbeshy[19]  analysed  the  problem  of  heat  transfer  over  an  exponentially  stretching  sheet  with  suction  confined  only  to 
viscous  fluid.  Subsequently,  Sanjayanand  and  Khan  [19,  20]  extended  the  work  of  Elbashbeshy  [18]  to  visco  -  elastic  fluid 
flow,  heat  and  mass  transfer  over  an  exponentially  stretching  sheet. 

The  investigation  of  exact  solutions  to  the  problem  of  nonlinear  equations  plays  an  important  role  in  the  study  of 
the  nonlinear  physical  phenomena.  Firstly,  the  Homotopy  Analysis  Method  [21,  22]  was  proposed  by  Liao  in  1992. 
Further,  the  HAM  was  developed  and  improved  by  Liao  for  nonlinear  problems  in  [23].  In  this  process,  Sajid  and  Hayet 
[24]  discussed  the  influence  of  thermal  radiation  on  the  boundary  layer  flow  due  to  an  exponential  stretching  sheet. 
Recently,  Hymavathi.  T  [25]  et.  al  analysed  the  flow  and  heat  transfer  in  visco-elastic  fluid  through  a  porous  medium  over 
an  exponentially  stretching  sheet.  They  used  Homotopy  Analysis  Method  (HAM)  to  solve  the  problem  analytically. 

In  view  of  the  importance  of  a  porous  medium,  we  investigate  to  study  of  the  heat  transfer  of  visco-elastic  fluid 
flow  over  exponentially  stretching  sheet  in  porous  medium  with  viscous  dissipation.  Here  we  employed  Homotopy 
Analysis  Method  to  solve  converted  non-linear  differential  equations  and  discuss  the  effect  of  visco  -  elastic  parameter 
(ki),  porosity  parameter  (k2),  Prandtl  number  (Pr)  and  Eckert  number  (Ec)  on  the  velocity  and  temperature  profiles. 
Comparison  of  the  present  analysis  is  also  made  with  the  existing  results  in  the  literature  and  is  seen  in  good  agreement. 

MATHEMATICAL  FORMULATION 


We  consider  steady  state  two-dimensional  incompressible  visco-elastic  fluid  over  exponential  stretching  sheet 
through  porous  medium  with  viscous  dissipation.  The  flow  will  be  generated  by  applying  two  equal  and  opposite  force 
along  x-axis  such  that  the  origin  is  fixed  and  the  flow  confined  to  y>0.  The  sheet  is  stretched  with  speed  varying 
exponentially  with  a  distance  from  the  slit.  We  take  x-axis  is  along  the  surface  and  y-  axis  being  perpendicular  to  it. 
The  governing  equations  for  the  flow  are  defined  as: 


du  dv 

- 1 - 

dx  dy 
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Where  9  is  the  kinematic  viscosity,  k0  is  the  elastic  parameter  and  k  is  the  permeability  of  porous  medium. 
Where  k  is  the  thermal  conductivity,  p  is  the  density,  T  is  the  temperature,  Cp  is  the  specific  heat  at  constant  pressure  and 
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is  viscous  dissipation. 


The  boundary  conditions  for  the  velocity  fields  are  of  the  form: 
u  =  UJx)  =  U0  e'1,  v  =  0,  T  =TW=  Tx  +  T0  e'1  city  =  0 
u  =  0,  uy  =  0,  T  =  Too,  as  y  — >  oo 


(4) 


where  Tw  is  the  temperature  of  the  sheet,  T-,  is  the  temperature  of  the  fluid  far  away  from  the  sheet  and  Z  is  the 
characteristic  length. 

SOLUTION  OF  THE  PROBLEM 

The  velocity  components  u  and  v  in  terms  of  stream  function  \|/(x,  y)  can  be  written  as: 


By/  By/ 

u  =  — — ,  v  = - 

By  Bx 


For  solving  momentum  equation,  introduce  a  similarity  variable  q  such  that 


11  =  y. 


I U 


2vl 


0  x/21 
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(6) 


y/(x, y)  =  ^2vlU0f(x,i])ex/21 , 0(tj) 


T  -Trx 
T-T 


(7) 


Here /is  the  dimensionless  stream  function  and  considering  if  (x,  tj)  =  f(rj).  Using  (5)  to  (7),  the  equations  (2)  & 
(3)  results  in  a  fourth  order  non-  linear  ordinary  differential  equation  of  the  form 


2fa-rr=r'-k1 


-k2f' 


(8) 


and 


O' Pr  (2  f'O-  fO')  +  k2  Pr  Ecf'2  +  Pr  Ecf"2  =  0 


(9) 
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,  Ku 

Where  k,  = - 

1  vl 


Pr 


Ec 


Uou, 
Ac  „ 


are  the  visco-elastic  parameter,  porosity 


parameter,  Pr  and  1 1  number  and  Eckert  number  respectively. 
Subject  to  the  boundary  conditions  are 


/  =  0 ,  /'=  1,  #  =  1 ,  at  q  =  0 

/'=  f"=0  =  0  asp— >oo  (10) 


ANALYTICAL  SOLUTION  (HAM) 

In  this  section,  we  employ  HAM  to  solve  the  equations  (8)  to  (9)  subject  to  the  boundary  conditions  (10). 
We  choose  the  initial  guesses  f0,  do  of  /,  9  in  the  following  form: 


fM= l-e"1  (ID 

OM^e-71  (12) 

The  linear  operators  are  selected  as 

£,(/)=/"'-/’  (13) 

L2{6)=0"-0  (14) 

which  have  the  following  properties 

L1(ci  +  C2e_'?+C3e'?)=  0  (15) 

L2(c4e"  +  C5e'1)=0  (16) 


Where  C  /  i  =  1  to  5  )  are  arbitrary  constants.  If  q  G  [() ,  1  ]  is  the  embedding  parameter,  ti  x  and  fl ,  are  the  non  - 
zero  auxiliary  parameters  and  H  j  (77),  H ,  (/))  are  the  auxiliary  functions,  then  we  can  construct  the  zeroth  -  order 
deformation  equations 

(1  - q)u q)-fo (17))  =qfrl  H1  {t])N]  (/ (77, q ))  ( 17) 

(l  -  q )L2  (< 9(jj , q)-0 0 (77))  =qh2  H2 {jj)N2  (/(t 7, q), q))  (18) 


subject  to  the  boundary  conditions 

f(0,qhh  f'{0,q)=l,  f'{oo,q)=0 

0{0,q)=l  0'{°°,q)=O 

Where 
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f(r/,l)=  f(i/) 
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(20) 

(21) 

(22) 


Thus,  as  q  increases  from  0  to  1,  f(r/,q )  varies  from  f{l/,0)  to  f{r/)  and  0{f/,q)  varies  from  0(77, 0)to 
0(77)  .  Then  expanding  f  (77 ,q\d{r/ ,q)  using  Taylor’s  theorem  with  respect  to  q,  we  obtain 

f(^q)=fo(7l)+YJfMqm  (23) 

777=1 


0(0;  4) = £0(0)+ IX  (0)4 m 

772=1 


Where 
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The  auxiliary  parameters  are  properly  chosen  so  that  series  (23)  and  (24)  converges  at  q  =1  and  thus 


(24) 


(25) 


(26) 


/(0)=/o(0)+Z/,„(0) 

777=1 


(27) 


0(0) = #0  (0) + X  (0) 

777=1 

(28) 

The  resulting  problems  at  the  mth  order  deformation  are 

A  La, (0)] = (0) . 

(29) 

A  k,  (0 )  -  xA-1  (0 )]  =  KK  (0 ) 

(30) 

Subject  to  the  boundary  conditions 

L  (0) = 0,  /■„  (0) = 0,  /■„  H = 0,  e,  (0) = a.e m  H = 0 

(31) 
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k= 0 


k= 0 


m- 1  m— 1  A  /?? — 1  1  m-\  Q  m-1  ' 

V  A-=0  ^k=G  ^k= 0  7 


(32) 


^  1  ,  2 
<  (v)  =  Cl  +  Pr  X  /m-1-/  /  +  £c  Pr  /m-1  +  Pr  ^2/„-l 


(=0 


(33) 


and 


f 0, 777  <  1, 

Xm~  \  *  , 

[  1,  m  >  1 

We  choose  the  auxiliary  function  as  follows: 


(34) 


//,(/7)=l  ,tf2(/7)=l.  (35) 

If  we  let  /„;(// M„,  (//)  as  the  special  solutions  of  the  m,h  order  deformation  equations,  the  general  solutions  are 

given  by 


fmi’l)=fM+ C1  +  C2e-'?  +  C3e'? 

^(7)=<?:(7)+C4e-’  +  C5e’ 

Where  the  integral  constants  C,  (/  =  1  to  5)  are  determined  using  the  boundary  conditions  (31). 


(36) 

(37) 


The  linear  non-homogenous  equations  (29)  and  (30)  are  solved  using  MATHEMATICA  software  one  after  the 
other  by  considering  m  =  1,2, . 

CONVERGENCE  OF  HAM  SOLUTION 

Liao  [16]  showed  that  for  an  analytic  solution  obtained  by  HAM,  its  convergence  and  rate  approximation  strongly 
depending  upon  the  auxiliary  parameters  ft  t  and  ft  ,  .  If  these  parameters  are  chosen  properly,  then  the  solution  is 

effective.  Hence,  fl  -curves  are  plotted  at  25th  order  approximation  in  order  to  obtain  the  suitable  ranges  for  -1.15  <  ft  {  <  - 
0.1,  -1.1  <  ftl  2—  -0.09,  with  convergence  of  li  value  is  h  =  -0.8.  Table  1  shows  the  convergence  of  the  solutions  with 
increasing  order  of  approximations. 

RESULTS  AND  DISCUSSIONS 

A  comprehensive  parametric  study  is  considered  for  flow  and  heat  transfer  of  visco-elastic  fluid  through  porous 
medium  with  viscous  dissipation  over  exponentially  stretching  sheet.  The  approximate  analytical  solutions  are  obtained  by 
the  Homotopy  Analysis  Method  for  physical  parameters  like  visco-elastic  parameter,  porous  parameter,  Pr  and  tl  number 
and  Eckert  number.  A  systematic  study  is  performed  to  analyze  the  effect  of  various  values  for  physical  parameters 
visco-elastic,  porosity,  pr  and  tl  and  Eckert  number  on  velocity  and  temperature  through  graphically. 
We  now  discuss  the  variations  of  the  physical  quantities  of  engineering  importance  that  is  the  local  skin  friction  coefficient 
-f  "(0)  and  local  Nusselt  number  -0’(0).  The  convergence  values  of  - f”(0 )  and  -8’(0)  displayed  in  Table  1. 
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Table  1:  Convergence  of  HAM  Solution  for  different 
Orders  of  Approximations  when 
kl  -  0.1,  k2  =0.1,  Ec  =  0.1,  Pr  =  1.0 


Order 

-/"( 0) 

-e\  0) 

5 

1.424645 

1.203512 

10 

1.427125 

1.200031 

15 

1.427156 

1.199619 

20 

1.427156 

1.199521 

25 

1.427156 

1.199492 

30 

1.427156 

1.199483 

35 

1.427156 

1.199481 

40 

1.427156 

1.199481 

Table  2:  Comparison  of  -0’(O)  for  different  Values 
of  k!  with  k2  =  0.5,  Pr  =  1.0,  Ec  =  0.25 


K, 

Mahentesh  M 

Present  Study 

0.0 

1.12221 

1.122209 

0.1 

1.09631 

1.0963096 

0.2 

1.06551 

1.065513 

From  the  figure  1  shows  the  admissible  values  of/”  (0)  and  9’  (0)  at  25th  order  approximation  on  velocity  and 
temperature  profiles. 

From  the  figures  2  and  3,  shows  the  effect  of  visco-elastic  parameter  on  velocity  and  temperature. 
It  can  be  found  that  the  increase  in  visco-elastic  parameter  decreases  the  velocity  as  well  as  increase  in  temperature. 
This  is  due  to  the  fact  that  elastic  property  in  visco-elastic  fluid  reduces  the  frictional  force. 


Figure  1:  25th  Order  Approximations  for  hi,  h2  Curves 
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k2=0.1  R=1.  Bc=0.1 


1 

Figure  2:  The  Effect  of  Visco-Elastic  Parameter  on  Velocity 
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Figure  3:  The  Effect  of  Visco-Elastic  Parameter  on  Temperature  Distribution 

From  figures  4  and  5,  shows  the  effect  of  porosity  parameter  on  velocity  and  temperature  distribution. 
It  can  be  found  that  the  increase  in  porosity  parameter  decreases  the  velocity  as  well  as  an  increase  in  temperature. 
This  is  due  to  increase  in  porosity  parameter  leads  to  the  enhanced  deceleration  of  the  flow,  thickening  of  the  thermal 
boundary  layer  occurs  due  to  the  increase  of  non-Newtonian  visco-elastic  normal  stress  and  the  effect  of  dimensionless 
porous  parameter  becomes  smaller  as  porous  parameter  increases. 


k1=  0.1  .  Pr=  1.  Ec=  o.i 


i 

Figure  4:  The  Effect  of  Porosity  Parameter  on  Velocity 
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kl  =  0.1  .  Pr=  1.  Ec=  0.1 


*1 

Figure  5:  The  Effect  of  Porosity  Parameter  on  Temperature 

From  the  figures  6  and  7,  shows  the  effect  of  Pr  and  tl  number  and  Eckert  number  on  temperature  distribution. 

It  can  be  found  that  increase  in  Pr  and  tl  number  decreases  the  temperature  distribution  and  increases  in 

Eckert  number  increases  the  temperature  distribution.  This  is  because  there  would  be  a  decrease  of  the  thermal 
boundary  layer  thickness  with  the  increase  of  values  of  Pr  and  tl  number  Pr.  The  increase  of  Pr  and  tl  number  means  slow 
rate  of  thermal  diffusion. 


0  2  4  6  8  10 

*1 


Figure  6:  The  Effect  of  Pr  and  tl  Number  on  Temperature  Distribution 
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Figure  7:  The  Effect  of  Eckert  Number  on  Temperature 
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CONCLUSIONS 

In  this  study,  Homotopy  Analysis  Method  is  employed  to  analyse  the  flow  and  heat  transfer  characteristic  of 
visco-elastic  fluid  over  exponentially  stretching  sheet  in  a  porous  medium  with  viscous  dissipation.  Analytical  solutions 
obtained  by  the  Homotopy  Analysis  method  (HAM)  are  seen  in  good  agreement  with  existing  results  in  the  literature. 
The  effect  of  various  parameters  on  velocity  and  temperature  is  analyzed  through  graphs.  The  results  are  summarized  as 
follows. 

•  The  effect  of  visco-elastic  and  porosity  parameter  on  horizontal  velocity  are  the  same. 

•  The  effect  of  Pr  and  tl  number  is  to  decrease  the  thermal  boundary  layer  thickness. 

•  The  effect  of  Eckert  number  increases  the  temperature  distribution 

•  The  visco-elastic  liquids  having  low  viscous  dissipation  must  be  chosen  for  effective  cooling  of  stretching  sheet. 
REFERENCES 

1.  Sakiadis,  B.  C.,  Boundary  Layer  Behaviour  on  Continuous  Solid  Surface:  1  -  Boundary  Layer  Equations  for  Two  Dimensional 
and  Axisymmetric  Flows,  AIChE.  J.,  7  (1961 ),  1,  pp.  26-28. 

2.  L.  J.  Crane,  Flow  past  a  stretching  plate,  Z.  Angew.  Math.  Phys.  21  (1970)  645-647. 

3.  V.  Kumaran,  G.  Ramanaiah,  A  note  on  the  flow  over  a  stretching  sheet,  Acta  Mech.  116  (1996)229-233. 

4.  K.  Vajravelu,  J.  R.  Cannon,  Fluid  flow  over  a  nonlinear  stretching  sheet,  Appl.  Math.  Comput.  181  (2006)  609-618. 

5.  H.  1.  Andersson,  MHDflow  of  a  viscoelastic  fluid  past  a  stretching  suiface,  Acta  Mech.  95  (1992)  227-230. 

6.  T.  Hymavathi,  B.  Sankar,  Heat  transfer  in  MHD  visco-elastic  fluid  flow  over  non-isothermal  stretching  sheet,  Journal  of  Appl. 
Mathematics  and  Computation  215(2009)2045-2054. 

7.  Cheng  P,  Minkowycz  WJ  (1977)  Free  convection  about  a  vertical  plate  embedded  in  a  porous  medium  with  application  to  heat 
transfer  from  a  dick.  J.  Geo.  phys.  Res  28:2040-2048 

8.  Subhas,  A.,  Veena,  P.,  Viscoelastic  Fluid  Flow  and  Heat  Transfer  in  a  Porous  Medium  over  a  Stretching  Sheet,  Int.  J. 
Nonlinear  Mech.  33  (1998),  3,  pp.  531-540. 

9.  Prasad,  K.  V.,  Abel,  M.  S.,  and  Datti,  P.  S.  (2003).  Diffusion  of  chemically  reactive  species  of  a  non-Newtonian  fluid  immersed 
in  a  porous  medium  over  a  stretching  sheet,  Int.  J.  Non-  Linear  Mech.,  38,  651-657. 

10.  Krishnambal,  S.,  and  Anuradha,  P.  (2006).  Effect  of  radiation  on  the  flow  of  a  visco-elastic  fluid  and  heat  transfer  in  a  porous 
medium  over  a  stretching  sheet,  J.  Appl.  Sci.,  6(14),  2901-2906. 

11.  M.  Subhas,  Abel,  S.  K.  Khan  and  K.  v.  Prasad,  Convective  heat  and  mass  transfer  in  a  visco-  elasticfluid  flow  through  a 
porous  medium  over  a  stretching  sheet,  Int.  Journal  of  Numerical  Methods  for  Heat  &  Fluid  flow,  Vol.ll  No.8,2001,  pp.779- 
792. 

12.  K.  A.  Yili,  Viscous  and  Joule  heating  effects  on  Non-Darcy  MHD  natural  convection  flow  overa  permeable  sphere  in  porous 
media  with  internal  heat  generation,  Int.  Commun.  Heat  Mass, 27(4),  pp.  591-600,  2000. 

13.  R.  M.  South,  S.  K.  Khan,  M.  S.  Abel,  K.  V.  Prasad,  Heat  and  mass  transfer  in  a  visco-elastic  fluid  over  an  accelerating  surface 
with  heat  source/sink  and  viscous  dissipation,  Heat  Mass  Transfer,  38,  pp.  213-220,  2002. 


Impact  Factor  (JCC):  6.2037 


NAAS  Rating:  3.76 


A  Homotopy  Analysis  Method  to  Flow  and  Heat  Transfer  Characterestics  11 

of  a  Visco-Elastic  Fluid  through  Porous  Medium  over  an 
Exponentially  Stretching  Sheet  with  Viscous  Dissipation 

14.  C.  H.  Chen,  Combined  heat  and  mass  transfer  in  MHD  free  convection  from  a  vertical  surface  with  ohmic  heating  and  viscous 
dissipation,  Int.  J.  Eng.  Sci.,  42,  pp.  699-713,  2004. 

15.  Mahantesh  M.  Nandeppanavar,  M.  Subhas  Abel,  K.  Vajravelu,  Flow  and  heat  transfer  characterestics  of  a  visco-  elastic  fluid 
in  a  porous  medium  over  an  impermeable  stretching  Sheet  with  viscous  dissipation,  Int.  J.  of  Heat  and  Mass  transfer 
53(2010)7407-4713. 

16.  W.  C.  Troy,  E.  A.  Overmanii  and  G.  B.  Ermentrout,  Uniqueness  of  flow  of  a  second  order  fluid  past  a  stretching  sheet, 
Quarterly  of  Applied  Mathematics,  Volume  XLIV,  Number4,  January  1987. 

17.  Elbashbeshy,  E.  M.  A.,  Heat  Transfer  Over  an  Exponentially  Stretching  Continuous  Surface  with  Suction,  Arch.  Mech.,  53 
(2001),  6,  pp.  643-651. 

18.  Rajalaxmi,  V,  Mohankumar,  G.,  &  Ramanathan,  K.  Effectiveness  Of  Plantar  Fascia  Stretching  Ks  Contrast  Bath  Combined 
With  Ultrasound  In  Plantar  Fasciitis. 

19.  Khan,  S.  K.,  Sanjayanand,  E.,  Viscoelastic  Boundary  Layer  Flow  and  Heat  Transfer  over  an  Exponentially  Stretching  Sheet, 
Int.  J.  Heat  Mass  Transfer,  48  (2005),  8,  pp.  1534-1542. 

20.  Sanjayanand,  E.,  Khan,  S.  K.,  On  Heat  and  Mass  Transfer  in  a  Viscoelastic  Boundary  Layer  Flow  over  an  Exponentially 
Stretching  Sheet,  Int.  J.  Therm.  Sci.,  45  (2006),  8,  pp.  819-828. 

21.  Liao  SJ.  The  proposed  homotopy  analysis  technique  for  the  solution  of  nonlinear  problems  (Ph.  D.  thesis),  Shanghai  Jiao  Tong 
University;  1992. 

22.  S.  J.  Liao,  A  new  branch  of  solutions  of  boundary-layer  flows  over  an  impermeable  stretched  plate,  International  Journal  of 
Heat  &  Mass  Transfer  48(2005),  2529  -  2539. 

23.  Liao  SJ.  Beyond  perturbation:  introduction  to  the  homotopy  analysis  method.  Boca  Raton:  Chapman  and  Hall/CRC  Press; 
2003. 

24.  Wu  W,  Liao  SJ.  Solving  solitaiy  waves  with  discontinuity  by  means  of  the  homotopy  analysis  method.  Chaos  Solitons  Fractals 
2005;26:177-85. 

25.  Sajid,  M..  Hayat,  T.,  Influence  of  Thermal  Radiation  on  the  Boundary  Layer  Flow  Due  to  an  Exponentially  Stretching  Sheet, 
Int.  Comm.  Heat  Mass  Transfer,  35  (2008),  3,  pp.  347-356. 

26.  T.  Hymavathi  and  B.  Akkayya,  Homotopy  Analysis  method  to  flow  and  heat  transfer  flow  in  visco-elastic  fluid  flow  in  porous 
medium  over  exponentially  stretching  sheet,  lilt.  Journal  of  Computational  and  Appl.  Mathematics,  Vol.l  lNo.2(2016),  pp.119- 
128. 


www.tjprc.ors 


editor@tjprc.  org 


